> 
o 

00 



o 



- 1—^ 

X 



Collimated Energy- Momentum 
Extraction from Rotating Black Holes 
in Quasars and Microquasars Using the 
Penrose Mechanism 



Reva Kay Williams* t 1 

(N 



' University of Florida, Gainesville, FL 32611 
^Bennett College, Greensboro, NC 27401 



Abstract. For almost four decades, since the discovery of quasars, mounting observa- 
tional evidence has accumulated that black holes indeed exist in nature. In this paper, 
I present a theoretical and numerical (Monte Carlo) fully relativistic 4-D analysis of 
Penrose scattering processes (Compton and 77 — > e~e + ) in the ergosphere of a super- 
massive Kerr (rotating) black hole. These model calculations surprisingly reveal that 
the observed high energies and luminosities of quasars and other active galactic nuclei 
(AGNs), the collimated jets about the polar axis, and the asymmetrical jets (which 
' can be enhanced by relativistic Doppler beaming effects) all are inherent properties 

of rotating black holes. From this analysis, it is shown that the Penrose scattered 
escaping relativistic particles exhibit tightly wound coil-like cone distributions (highly 
collimated vortical jet distributions) about the polar axis, with helical polar angles of 
escape varying from 0.5° to 30° for the highest energy particles. It is also shown that 
the gravitomagnetic (GM) field, which causes the dragging of inertial frames, exerts 
a force acting on the momentum vectors of the incident and scattered particles, caus- 
ing the particle emission to be asymmetrical above and below the equatorial plane, 
jJJ ■ thus appearing to break the equatorial reflection symmetry of the Kerr metric. When 

the accretion disk is assumed to be a two-temperature bistable thin disk/ion corona 
(or torus = advection-dominated accretion flow), energies as high as 54 GeV can be 
attained by these Penrose processes alone; and when relativistic beaming is included, 
energies in the TeV range can be achieved, agreeing with observations of some BL Lac 
objects. When this model is applied specifically to quasars 3C 279 and 3C 273, their 
observed high energy luminosity spectra can be duplicated and explained. Moreover, 
this energy extraction model can be applied to any size black hole, irrespective of 
the mass, and therefore applies to microquasars as well. When applied specifically to 
microquasar GRS 1915+105 the results are consistent with observations. 
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I INTRODUCTION 



Astrophysical jets are one of the most poorly understood phenomena today. It 
is clear that they are present where gravitational accretion or contraction and 
magnetic fields exist. We observe these jets in quasars and microquasars due to 
supermassive and stellar size black holes, respectively. They are also present in 
contracting protostars. Perhaps, understanding the mechanism responsible for jets 
appearing from the energy source of black holes, where the gravitational field is 
dominant, we can understand their appearences associtated with photostars. At 
present there are two popular trains of thought associated with jets in black holes: 
one is that the jets are inherent properties of geodesic trajectories in the Kerr [1] 
metric of a rotating black hole, and thus, can be described by Einstein's general 
theory of relativity; and the other is that the accretion disk and its magnetic 
field through magnetohydrodynamics (MHD) are producing the jets. Perhaps it 
could be a combination of the two, with gravity controlling the flow near the event 
horizon, and MHD controlling the flow at distances farther away. In this paper, 
an anaylsis of the Penrose mechanism [2] is presented to describe gravitational- 
particle interactions close to the event horizon in the subparsec regime. In this 
fully general relativistic description, the jets are produced and controlled by gravity 
alone, without the necessity of an external magnetic field. 

II MODEL 

The model consists of a supermassive 1O 8 M Kerr (rotating) black hole plus 
particles from an assumed relativistic bistable thin disk/ion corona [or torus = 
advection-dominated accretion flow (ADAF)]: two-temperature [separate tempera- 
tures for protons (~ 10 12 K) and electrons (~ 10 9 K)] accretion flow. The Penrose 
mechanism is used to extract rotational energy-momentum by scattering processes 
inside the ergosphere (r erg ~ 2M, where c = G — 1). See Williams [3] for a de- 
tailed description of the model. The "quasi-Penrose" processes investigated are 
(a) Penrose Compton scattering (PCS) of equatorial low energy radially infalling 
photons by equatorially confined (Q e = 0) and nonequatorially confined (Q e ^ 0) 
orbiting target electrons, at radii between the marginally bound (r mb ~ 1.089M) 
and marginally stable (r ms ~ 1.2M) orbits; (b) Penrose pair production (PPP) 
(iP — > e~e + p) at r m t,; (c) PPP (77 — > e~e + ) by equatorial low energy radi- 
ally infalling photons and high energy blueshifted (by factor ~ 52) nonequatorially 
confined 7-rays at the photon orbit (r ph ~ 1.074M). Note, the target particles are 
initially in bound (marginally stable or unstable) trapped orbits, trapped in the 
sense of possibly having no other way of escaping save these Penrose processes [4,3]. 
Note also that, as the nonequatorially confined target particle passes through the 



equatorial plane, Q, a constant of motion as measured by an observer at infinity 
[5,3], equals P@, where P® is the polar coordinate momentum of the particle. 

Ill METHOD 

Monte Carlo computer simulations of up to ~ 50, 000 scattering events of infalling 
accretion disk photons (normalized to a power-law distribution) are executed for 
each quasar. Energy and momentum (i.e., 4- momentum) spectra of escaping par- 
ticles (7's, e~e + 's), as measured by an observer at infinity, are obtained (where 
7 = a general photon, i.e., of any energy value). The following ingredients are 
used: (1) General relativity [the Kerr metric spacetime geometry yields equatori- 
ally and nonequatorially confined ("spherical-like" [6]) particle orbits and escape 
conditions, conserved energy and angular momentum parameters, and transfor- 
mations from the Boyer-Lindquist coordinate frame (BLF) to the local nonrotat- 
ing frame (LNRF)]. Note, BLF is the observer at infinity [7]; LNRF is the local 
Minkowski (flat) spacetime. (2) Special relativity [in the LNRF, physical processes 
(i.e., the scatterings) are done; Lorentz transformations between inertial frames are 
performed; and Lorentz invariant laws are applied]. (3) Cross sections [application 
of the Monte Carlo method to the cross sections, in the electron rest frame for PCS, 
in the proton rest frame for PPP(7£> — > e~e + p), and in the center of momentum 
frame for PPP(77 — > e~e + ), give the distributions of scattering angles and final 
energies] . 

IV OVERALL RESULTS 

The energies attained are the following: (1) PCS: For the input (photon) energy 
range 5 eV to 1 MeV, the corresponding output energy range is ~ 15 keV to 14 
MeV. (2) PPP (7^ — ► e~e + p): No escaping pairs for radially infalling 7-rays (~40 
MeV), and no energy boost: implying that the assumption (negligible recoil energy 
given to the proton) made in the conventional cross section and, perhaps, the 
geometry of the scattering must be modified. It had been predicted [8] that pairs 
with energies (~1 GeV) can escape. (3) PPP (77 — ► e~e + ): For input (photon) 
energy range ~ 3.5 keV to 100 MeV, yields output (e~e + ) energy range ~ 2 MeV 
to 10 GeV (for BB), and higher up to ~ 54 GeV (for PL, input 7 ~ 2 GeV), where 
BB, PL = blackbody and power-law distributions, respectively, for the accretion 
disk protons that yield the neutral pion decays it — > 77 [9,10] to populate the 
photon orbit. 

The luminosity spectrum due to Penrose processes for the specific case of quasar 
3C 273 is plotted in Figure 1(a), along with the observed spectrum for comparison. 
The outgoing (escaping) luminosity spectrum produced by the Penrose scattered 
particles is given by [11] 

L C J C « 4ird 2 F° sc (erg/sHz) 



^And 2 hu csc f 1 f 2 ---f n (N^-Nr), 

where d is the cosmological distance of the black hole source; _F° SC is the flux of es- 
caping photons; iV* n and A^ ap are the emittance of incoming and captured photons, 
respectively; f n defines the total fraction of the particles that undergoes scattering 
[n = 2 for PCS and n = 5 for PPP (77 — ► e"e+)]. The values of fi,...,f n 
are the fitting factors, which make the Penrose calculated luminosities agree with 
observations for the specific case of 3C 273. Note, in the model calculations, if 
we let every particle scatter, and set /1 = fa ~ 1CT 2 , defining the fraction of the 
disk luminosity intersecting the scattering regime, with the remaining /„'s equal 
1, the continuum emissions (the top curves on Figures 1(a), labeled with numbers 
for specific cases of target and incident particles) are obtained; see [11] for details. 
The spectrum resulting from the PPP (77 — > e~e + ) is produced by letting the 
escaping pairs undergo "secondary Penrose Compton scattering" with low energy 
radially infalling accretion disk photons (= / 3 ). 

Thus, as one can see from Figure 1(a), the Penrose mechanism can generate the 
necessary luminosity observed, and the three model calculated regions of emission 
[due to PCS by equatorial confined targets (curve passing through nos. 1 — 7), 
by nonequatorially confined targets that cross the equatorial plane (curve passing 
through nos. 8 — 13), and PPP (77 — > e~e + ) (curve passing nos. 14 — 25)] are 
consistent with the three major regions of emission in all quasars and AGNs. Note, 
the target photons at the photon orbit can only exist in nonequatorially confined 
orbits [3]. 

Moreover, the observed spectra of microquasars (or galactic black holes), in gen- 
eral, appear not to have PCS emission by the nonequatorially confined target elec- 
trons, as well as the highest energy 7-ray emission due to PPP (77 — > e~e + ), 
indicating that these sources may not have an ion corona nor ADAF, which would 
be need to populate the orbits to generate such emission. 

V THE GRAVITOMAGNETIC FIELD 

The gravitomagnetic (GM) force field is the gravitational analog of a magnetic 
field. It is the additional gravitational force that a rotating mass produces on a test 
particle. The GM force is produced by the gradient of [3 GM = —uje$>, where uj is 
the frame dragging velocity and (3 GM is gravitomagnetic potential [12]. Analysis of 
the equations governing the trajectories of the particles shows that the GM force, 
which acts proportional to the momentum of a particle, alters the incoming and 
outgoing momentum parameters of the incident and scattered particles, resulting in 
asymmetrical distributions, thus, appearing to break the reflection symmetry of the 
Kerr metric above and below the equatorial plane [3,13,14]. Effects of the GM force 
acting on the PPP (77 — ► e~e + ) process can be discerned from Figures 1(d), 1(e), 
and 2. Note, only the distribution for (P^e > is shown in Figure 1. When half of 
the 2000 target photons are allowed to have (P 7 2)e > and the other half (P 7 2)e < 



of equal absolute values, the e~e + "jet/counter-jet" achieve a maximum ratio 
~ 3 : 1, favoring (P T )e > [13], as seen in Figure 2. Polar coordinate momentum 
distributions, (Pp h )e, for PCS escaping photons are displayed in Figure 3. Notice 
the effects of the GM force field causing the photon jet /counter-jet to vary from 
nearly symmetric to asymmetric for the different cases shown. Of these cases the 
largest ratio achieved is ~ 5 : 1 [Fig. 3(c)]. The direct cause of the asymmetry in the 
polar direction appears to be due to the severe frame dragging in the ergosphere in 
which the GM field lines are dragged in the direction that the black hole is rotating 
[13]. 

VI THE VORTICAL ORBITS PRODUCED 

It is found that the Penrose scattered particles escape along vortical trajectories 
collimated about the polar axis [3,11,16]. These distributions are fluxes of coil- 
like trajectories of relativistic jet particles escaping concentric the polar axis. The 
highest energy escaping particles have the largest P$ values, the near largest P® 
values, and smallest \P r \ values. Note, P r is negative (inward toward the polar axis) 
for most of the PCS photons, and positive for the e~e + pairs. The helical angles of 
escape (<5j) csc = 1 90° — ^ | of particle type i, for the highest energy scattered particles 
are (<5 ph )csc ~ 1° - 30° for PCS and (5 T ) CSC ~ 25° - 0.5° for the e"e+ pairs. The 
above characteristics of the escaping particles imply strong collimation about the 
polar axis, giving rise to relativistic jets with particle velocities up to ~ c [eg. cf. 
Figs. 1(b), 1(c), 1(f), and 2]. Note, such vortical trajectories and collimation are 
consistent with the findings of de Felice et al. from spacetime geometrical studies 
of general particle geodesies in a Kerr metric [18-20]. 

VII CONCLUSIONS 

From this model to extract energy-momentum from a black hole we can conclude 
the following: PCS is an effective way to boost soft x-rays to hard x-rays and 7-rays 
up to ~ 14 MeV. PPP (77 — > e~e + ) is an effective way to produce relativistic 
e~e + pairs up to ~ 54 GeV: This is the probable mechanism producing the fluxes 
of relativistic pairs emerging from cores of AGNs. These Penrose processes can 
operate for any size rotating black hole, from quasars to microquasars [17]. Overall, 
the main features of quasars: (a) high energy particles (x-rays, e~e + pairs, 7-rays) 
coming from the central source; (b) large luminosities; (c) collimated jets; (d) one- 
sided (or uneven jets), can all be explained by these Penrose processes to extract 
energy from a black hole. 

Moreover, it is shown here that the geodesic treatment of individual particle 
processes close to the event horizon in the subparsec regime, as governed by the 
black hole, is sufficient to described the motion of the particles. This finding is 
consistent with MHD through the statement made by de Felice and Zanotti [20], 
that the behavior of such individual particles on geometry (or gravity) -induced 



collimated trajectories is also that of the bulk of fluid elements in the guiding 
center approximation. In light of this, MHD should be incorporated into these 
calculations to describe the flow of the Penrose escaping particles away from the 
black hole, i.e., to perhaps further collimate and accelerate these jet particles out 
to the observed kpc distances. 
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FIGURE 1. (a) Comparing the theoretical spectrum with observations for 3C 273. The cal- 
culated PCS and PPP (77 — ► e~e+) luminosity spectra are represented by the solid squares 
and large solid dots, respectively. The observed spectra is indicated by the solid line. The upper 
curves with the solid squares and solid dots superimposed on the dotted line and the dashed 
line, respectively, for PCS and PPP (77 — ► e~e+), are the spectra calculated from this model. 
Superimposed on the lower solid line of the observations are solid squares and solid dots that have 
been adjusted to agree with observations. These adjustments depend on the f„'s values (see text), 
(b) and (c) PPP (77 — ► e~e+) at r p h = 1.074M: scatter plots showing momentum components 
(each point represents a scattering event). The radial momenta (P^) r vs - E T , the azimuthal 
momenta (P^)* (= L T ) vs. E T ; for the infalling photons E 7 i = 0.03 MeV, and for the target 
photons E 7 2 = 3.893 GeV. (d) and (e) PPP (77 — ► e"e+): (P T ) vs. E T ; for E 7l = 0.03 MeV, 
£ 7 2 = 13.54 MeV, (P 72 )e = 0.393 Mm e ; and for E lX = 0.03 MeV, E l2 = 3.893 GeV, 
(^72)0 = 113 Mm e ; respectively, (f) The velocity distribution vs. 7 e (= E T /m e c 2 ) for the 
same case as (d) above. 
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FIGURE 2. PPP (77 — ► e~e + ) at r p h = 1.074M: scatter plots showing momentum compo- 
nents (each point represents a scattering event); for the infalling photons E 7 i = 0.03 MeV, and 
for the target photons E l2 = 3.893 GeV, (P 7 2)e = ±113 Mm e [similar to Fig. 1(e); see text]. 
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FIGURE 3. PCS: scatter plots showing polar coordinate space momenta: (-Pp h )e [= (Qph) 1 ^ 2 ] 
vs. £p h , of the escaping PCS photons after 2000 events (each point represents a scattering event), 
at r m b — 1.089M. The various cases are defined by the following parameters: E p h, initial photon 

1/2 

energy; E e , the target electron orbital energy; Q c , defining the corresponding polar coordinate 
momentum (P e )e of the target electron; iV os , number of photons escaping, (a) E p h = 3.5 keV, 
E e = 0.539 MeV, Ql /2 = 0, N cs = 1637. (b) E ph = 0.03 MeV, E e = 0.539 MeV, Ql /2 = 0, 
iV os = 1521. (c) E ph = 0.15 MeV, Ql /2 = 0, iV os = 1442. (d) E ph = 0.15 MeV, E e = 1.297 MeV, 
qV 2 = ±2.479 Mm e , N es — 1628. (Note, due to a minor oversight leading to improper treatment 
in the computer simulation of the arccosine term in eq. (3.39) of [3], correct Figs. 3(a) and 3(b) 
presented here replace Figs. 7(a) and 3(c), respectively, of [3]. ) 



